Vanishing of the Dissipationless Spin Hall Effect in a Diffusive Two-Dimensional 

Electron Gas with Spin-Orbit Coupling 
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We propose a nonequilibrium Green's function approach to study the spin-Hall effect in a two- 
dimensional electron system with both the Rashba and Dresselhaus spin-orbit couplings. By taking 
into account the long-range electron-impurity scattering, the derived kinetic equations are solved 
numerically. It is found the vanishing of the total zero-temperature dissipationless spin-Hall effect, 
contributing from the intrinsic and disorder-mediated processes. This result has been examined in 
the wide ranges of spin-orbit coupling constants and electron density. 

PACS numbers: 72.10.-d, 72.25.Dc, 73.50.Bk 
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In the growing field of spintronics, which aims 
to the manipulation of spin degrees of freedom in 
semiconductors Q, the issue of creating spin polariza- 
tion of carriers in nonmagnetic semiconductors with spin- 
orbit coupling by employing only the electric field has 
attracted a great deal of theoretical and experimental 
attention. More recently, the non-dissipative spin cur- 
rents with the spin-polarization perpendicular to the 
flow direction and to the applied electric field have 
been predicted in p-doped bulk semiconductor l2| and 
two-dimensional (2D) systems with the Rashba [3j and 
Dresselhaus spin-orbit couplings 0,01 . This dissipation- 
less spin-Hall effect essentially arises from the dc-field- 
induced transition between the hclicity bands: the dc 
field causes the variation of the electron momentum, 
but leaves its spin unchanged^]. In experiment, the 
spin-Hall effects have been observed in n-doped bulk 
semiconductors @ and two-dimensional hole gas@. 

The disorder effects on the spin-Hall current have 
been extensively investigated for the two-dimensional 
electron systems with Rashba spin-orbit coupling. In 
Ref. 13, llCl lll| , it has been argued that the disorder can 
reduce the spin-Hall effect and the ballistic value of spin- 
Hall conductivity reaches for the weak scattering limit. 
However, based on the Kubo formalism, Inoue et ol. \l2\ . 
Dimitrova[l3| and Chalaev and Loss0| have shown that 
the spin-Hall effect disappears even when the disorder is 
weak. This conclusion has been confirmed in Refs. 0] 
and by employing the spin-density matrix method 
and the Keldysh approach, respectively. Further, using 
the nonequilibrium Green's function approach, Liu and 
Lei have demonstrated that, even for the case of long- 
range disorders, the dissipationless spin-Hall effect also 
vanishes 6]. Obviously, the nonexistence of dissipation- 
less spin-Hall effect still holds for the Dresselhaus spin- 
orbit coupling, since it relates to the Rashba coupling 
through a global unitary transformation^. It is com- 
monly believed that this complete cancellation of spin- 
Hall effect is not a consequence of any symmetries |l 7j 
and relates to the isotropy of the dispersion in the helic- 
ity basis. 



It is very interesting to examine the spin-Hall effect for 
2D electron systems with both the Rashba and Dressel- 
haus spin-orbit interactions. In such systems, the spin- 
orbit-coupled bands become anisotropic, leading to many 
interesting phenomena, such as the spin splitting[18|, 
spin precession and relaxation^J HH EI E3- an d 
anisotropic transport [2^] , etc. In particular, when the 
Rashba coefficient is turned to be equal to the Dressel- 
haus one, such system can be used as the new type of spin 
field-effect transistor (SFET) operated in the nonballis- 
tic regime [il]. However, the spin-Hall effect in the 2D 
system with the Rashba and Dresselhaus spin-orbit cou- 
plings has been studied only in the ballistic regime 0,0]. 
It has been demonstrated that the spin-Hall conductivity 
has a universal value, independent on the strength of the 
coupling, but its sign is determined by the relative ratio 
of the two couplings. 

In this letter, we develop a nonequilibrium Green's 
function approach to the spin-Hall effect in 2D electron 
systems with both the Rashba and Dresselhaus spin-orbit 
couplings. The electron- impurity scattering is consid- 
ered in the self-consistent Born approximation and the 
obtained equations are solved numerically. We find that 
the dissipationless spin-Hall conductivity vanishes at zero 
temperature. This result contrasts with the general belief 
that the anisotropy can lead to the nonzero dissipation- 
less spin-Hall conductivity in 2D electron system with 
spin-orbit coupling [l2j. 

We consider an effective Hamiltonian for a two- 
dimensional electron of momentum p = (p x ,p y ) and ef- 
fective mass to 



Ho — ~ 1- Hr + Hi 

2m 



(1) 



with the terms of Rashba spin-orbit coupling Hr |24j , 

Hr = a(p y cr x - p x <7 y ), (2) 
and Dresselhaus spin-orbit interaction|25j. 



(3(p y (7y -p x a x ). 



(3) 
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Here, a = (a x ,a y , a z ) are the Pauli matrices, and, the a 
and /3 are the coupling constants. After taking the local 
unitary transformation 

U(P) = ^ ( _ e U e L ) ( 4 ) 

with the x p being 

X P = axg[ap y - /3p x - i(ap x - /3p y )], (5) 

we can diagonalize this Hamiltonian in the helicity basis 
H = diag(ei(p), e 2 (p))- The £ M (p) = 1,2) are of the 
forms 

= + (-ireflc, (6) 

where, e RD = y/(ap y - [3p x ) 2 + (ap x - f3p y ) 2 . 

When the system is driven by a dc field E applied along 
the x-direction, the i-direction-polarized spin current for 



electron with momentum p can be defined as j z = [a ? ;V + 
vc7j]/4 with the electron velocity v = — i[r, H] |26(. In 
the helicity basis, the total z-direction-spin current along 
the y axis, corresponding to the spin-Hall current, can be 
written as 

J v = E f^^P) + <° 21 (P)1 = E ^ Re <Mp), (7) 
p p 

where p^^(p) is the distribution function with the sym- 
metry relation puip) = P2i(p)II3- Further, the spin- 
Hall conductivity is given by <j s h — J y /E. 

To derive the kinetic equations for the 2D systems 
with the spin-orbit coupling, we follow the procedure de- 
scribed in Ref. 0] . We first carry out the Dyson equations 
for the real-space less Green's functions in the spin ba- 
sis. While the Fourier and the unitary transformation 
(0} are performed, these equations reduce to the kinetic 
equations in the helicity basis, 



ji A + ieE . v p | G< + ^ • V P X P [G< a x ] - e RD [G<,a z } = £ r G< + £<G Q - G'S< - G<£ a 
here, the self-energy for the electron-impurity interaction with scattering matrix V(p — k) reads 
57><(p,T, t) = l -n % J2 \V(p - k)| 2 {aiG r '< + a 2 a x G r ><a x - ia 3 [a x , G r ><}} , 

k 

I 



(8) 



(9) 



with the center-of-mass and difference times T and t re- 
lating to the two times of the Green's functions ti and t 2 : 
T = {t\ + t-z) /2, t = ti — t 2 . In these equations, for short- 
ness, the arguments (k, T, t) of Green's functions G r,< are 
dropped. a;(i = 1,2,3) are the factors associated with 
the directions of the momenta, a% = 1 + cos(x P — Xk)i 
a 2 = 1 — cos(x p — Xk), a 3 = sin(x p — Xk)- We can see that 
in comparison with the case of only the Rashba termQ, 
here, the factor a, is dependent on the momentum angle 
complicatedly through \p- This fact leads to the kinetic 
equations analytically unsolvable. 
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To further simplify our treatment, we adopt the gener- 
alized Kadanoff-Baym ansatz (GKBA), which expresses 
the causality of the time development of the two-time 
propagators using its equal-time value [2^, |2j|. At 
the same time, we neglect the collision broadening of 
the retarded and advanced Green's functions and re- 
tain the quan tities only in the lowest order of gradient 
expansion |2jj. These simplifications correspond to the 
Boltzmann limit and in result, our formalism becomes 
quasiclassical. The kinetic equation for the distribution 
functions can be followed 



where 



dp 
df 



ieE 

PKP, T) — Y ' VpXp a ^ + %SRD ^ p ' ° z 



dp 
dT 



scatt J— oo 



/ dt' [S>G< +G<E> -£<G> -G>S<] (T,t')(t' ,T). 

J —oo 



(10) 



(11) 



Further, we assume that the dc field is sufficiently weak and only the linear response in the steady state needs to 
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be considered. Otherwise, since we are just interested 
in the dissipationless spin-Hall effect, it is sufficient to 
treat the nondiagonal distribution functions in the lowest 
order of electron- impurity scattering, i.e. in the order of 
(rii)°. Under these considerations, the solutions of the 
kinetic equations (fTUf) can be separated into two parts. 
The first part of solution associates with the driving term 
■|eE-V p Xp[poj fj with po being the distribution function 
in equilibrium po(p) = diag[np(ei(p)), np(£2(p))], and is 
a nondiagonal matrix in the spin space. It can be written 
as 

p$(p) = P2i } (p) = eE 4g ^* P WMp)] - ^fMp)]} , 

(12) 

and results in the nonzero spin-Hall conductivity at zero 
temperature^, , 



r (i) = J _v_ 

r sH- E 



-e/87r a > (3 
a = (3 
e/8n a < (3 



(13) 



This part of spin-Hall conductivity corresponds to 
the intrinsic spin-Hall effect and agrees with the pre- 
vious studies in the ballistic regime 0, 0|. Essentially, 
this dissipationless spin-Hall effect stems from the dc- 
field-induced transitions between two spin-orbit-couplcd 
bands. When a dc field is applied to the 2D system with 
spin-orbit coupling, electron obtains an additional mo- 
mentum, but its spin remains unchanged. Consequently, 
in the helicity basis, the transition between the spin- 
orbit-coupled bands occurs. It is obvious that all elec- 
trons in 2D systems join in this process. 



The second part of solution relates to the transport 
process and the corresponding diagonal distribution func- 
tions should be determined from the coupled equations 



dn F (E) 



dE 



dpa 



= nJ2 \V(P k)| 2 {[pg)(p) - pW(k)] ai A^ + [pg](p) - p§(k)]a 2 A^}, 



(14) 



with /Li = 1,2, p, = 3 — p, and A M „ = <5(e A1 (p) — e„(k)). And then, the real part of nondiagonal distribution can be 
given by 



Rep£> = Rep%> = J E I^P " k)| 2 a 3 (-l)"{A w b< 2 »( P ) - p$(k)] - A MA [p$(p) _ p g>(k)]}. 



(15) 



k/x=l,2 



This part of spin-Hall effect depends on the electron- 
impurity scattering and entirely vanishes in clean sam- 
ples. However, it is of order of (n^) and the disorder 
only plays an intermediate role. This disorder-mediated 
spin-Hall effect contributes from the electrons near the 
Fermi surface and, physically, can be understood as the 
following way. The electrons, participating in the lon- 
gitudinal transport, also can be scattered by impurities. 
This process, which changes the electron momentum and 
conserves its spin, yields the transition between two spin- 
orbit-coupled bands in the helicity basis. 

Unfortunately, the above coupled equations for diag- 
onal distribution functions only can be resolved numer- 
ically. We perform a numerical evaluation for a GaAs- 
based heterojunction, where the electron-impurity scat- 
tering is assumed to be long-range The momentum 
integration is calculated by the Gauss-Legendre scheme. 
The linear system equations, arising from the discretiza- 
tion of the integral equations for diagonal distribution 
functions, are solved by using the singular value decom- 
position. Here, we concern on the zero-temperature spin- 
Hall effect. Therefore, only the distribution functions at 



the Fermi surface need to be carried out. 

We know that, when only presence of the Rashba term, 
the dependence of the nondiagonal distribution function 
on the momentum angle can be expressed through the 
sine function 6]. However, in our case, the deviation from 
this rule is expected due to the anisotropic dispersion re- 
lations. In Fig. 1, we plot the real part of nondiagonal 
distribution functions, which could be observed by the 
experimental tools, such as Raman spectroscopy, angle 
resolved photoelectron spectroscopy (ARPES) etc. It can 
be seen that when the f3 closes to the value of a, a peak 
always appears near the angle 7r/4 or 5ir/4. With the 
ratio of (3/a approaching to the unit, the peak becomes 
more pronounced. The amplitude of Repi2 for the same 
value of (3/a increases with the parameter a. It is obvi- 
ous that such structures near 7r/4 and 57r/4 relate to the 
zero of erd at a = /?, producing a jump in the electron 
velocity. 

Further, the spin-Hall current has been carried out by 
substituting the Repi2 into Eq. (jZJ). The total spin-Hall 
current, coming from the above two parts of the kinetic- 
equation solutions, is found to vanish. We have exam- 
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FIG. 1: The real part of total nondiagonal distribution func- 
tions pi2 near the angles 7r/4 and 57r/4. The electron density 
is 5 x 10 n cm- 2 . 



of spin- Hall current, associating with the transport, is 
expected to deviate from the linear rule. Hence, the non- 
dissipative spin-Hall effect can be found for the strength 
of dc field larger than about 1 kV/cm. 

For 2D electron systems with both the Rashba and 
Dresselhaus spin-orbit couplings, we even have calculated 
the longitudinal spin conductivity a% x , defined in Ref . . 
Its value also has been found to be zero. 

In summary, we construct a nonequilibrium Green's 
function approach to the spin-Hall effect of 2D electron 
systems with both the Rashba and Dresselhaus spin-orbit 
interaction. Taking into account the long-range electron- 
impurity scattering, we have resolved the kinetic equa- 
tions numerically. It is found that, although the dis- 
persion relation of such systems exhibits an anisotropic 
character, the total dissipationless spin-hall conductivity 
vanishes. 
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W. Xu. This work was supported by the National Sci- 
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State Basic Research Project, and the Youth Scientific 
Research Startup Foundation of SJTU. 



ined this result in the wide range of parameters a and 
(3 ~ 1 — 10 x lCP n eVm. The typical electron concen- 
trations are taken to be 1 — 10 x 10 n cm _2 . We have 
seen, although the nondiagonal distributions are depen- 
dent on the coupling constants and electron density, the 
total spin-Hall effect in the estimated errors always be- 
comes unobservable. 

This complete cancellation of dissipationless spin-Hall 
effect supports the phenomenological argument in Ref. 
|l4|: the vanishing of spin-Hall current comes from the 
fact that the uniform spin current is the total time deriva- 
tive of the magnetization. Here, when the presence of 
both the Rashba and Dresselhaus spin-orbit couplings, 
the corresponding magnetization operator can be taken 
as pS x + aS y with Si being the spin operator. We should 
note that the nonexistence of dissipationless spin-Hall ef- 
fect is not a consequence of any symmetries and is not 
necessarily suppressed to zero in general [13 • ^ depends 
on the dispersion relation and wave functions of the con- 
sidered systems. For example, we have carried out the 
spin-Hall conductivity for 2D hole systems and found a 
nonzero dissipationless spin-Hall conductivity. 

However, the dissipationless spin-Hall effect is ex- 
pected to be observed in the nonlinear regime. We can 
see from the kinetic equations, the intrinsic part of spin- 
Hall conductivity always linearly depends on the strength 
of dc fields. However, it is well known that the trans- 
port of 2D electron systems becomes nonlinear when the 
strength of dc field is of order of lkV/cm. At that 
time, the dc-field dependence of disorder-mediated part 
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